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Abstract: We present an extraction of the valence transversity parton distributions based 
\^ on an analysis of pion-pair production in deep-inelastic scattering off transversely polar- 

ized targets. Recently released data for proton and deuteron targets at HERMES and 

COMPASS permit a flavor separation of valence transversities. The present extraction 
,__( is performed in the framework of collinear factorization, where dihadron fragmentation 

functions are involved. The latter are taken from a previous analysis of electron-positron 
• . annihilation measurements. 
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1 Introduction 

The distribution of quarks and gluons inside liadrons can be described by means of parton 
distribution functions (PDFs). In a parton model picture, PDFs describe combinations of 
number densities of quarks and gluons in a fast-moving liadron. The knowledge of PDFs is 
crucial for our understanding of QCD and for the interpretation of high-energy experiments 
involving hadrons. At leading twist, the quark structure of spin-half hadrons is described 
by three PDFs: the unpolarized distribution the longitudinal polarization (helicity) 

distribution gi{x), and the transverse polarization (transversity) distribution hi{x) [1-4]. 
From the phenomenological point of view, the unpolarized PDFs are well-known, as can be 
evinced by the large number of parametrizations available (see, e.g., Ref. [5] and references 
therein). Apart from giving us invaluable information about the structure of nucleons, they 
have a fundamental importance for the interpretation of measurements in any hadronic 
colliders, e.g., the LHC. The helicity PDFs are known to some extent, see e.g. Refs. [6-9]. 
On the other hand, the transversity distribution is poorly known (see, e.g., Refs. [10] and 
references therein). This is mainly due to the fact that transversity can be measured only 
in processes with two hadrons in the initial state [11], e.g. proton-proton collision, or one 
hadron in the initial state and at least one hadron in the final state, e.g. semi-inclusive DIS 



Combining data from HERMES [12] and COMPASS [13] on polarized SIDIS with one 
hadron in the final state, together with data from Belle [14] on almost back-to-back emission 
of two hadrons in e^e~ annihilations, the transversity distribution was extracted for the first 
time by the Torino group [15]. The main difficulty of such analysis lies in the factorization 
framework used to interpret the data, since they involve Transverse Momentum Dependent 



(SIDIS). 
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PDFs (TMDs). In spite of exceptional progress in the understanding of TMDs [16-21], we 
have still limited information on their evolution equations, which are needed when analyzing 
measurements at very different scales. 

In this paper, we extract the transversity distribution for the valence combination of 
up and down quarks, applying for the first time an approach based on standard collinear 
factorization. We use data on SIDIS with two hadrons detected in the final state, where the 
transversity distribution is combined with the so-called Dihadron Fragmentation Functions 
(DiFFs) [22-24]. The collinear framework allows us to keep under control the evolution 
equations of DiFFs [25]. 

In Sec. 2, we summarize the theoretical framework for two-hadron SIDIS. The parametriza- 
tion of the valence transversity and its error analysis is described in Sec. 3. The results 
are discussed in Sec. 4. We finally illustrate the possible applications and extensions of our 
analysis and draw our conclusions in Sec. 5. 



2 Theoretical framework for two-hadron SIDIS 

We consider the process 

i{k) + N{P) ^ i{k') + Fi(Pi) + H2iP2) + X , (2.1) 

where i denotes the beam lepton, N the nucleon target. Hi and H2 the produced hadrons, 
and where four-momenta are given in parentheses. We work in the one-photon exchange 
approximation and neglect the lepton mass. We denote by M the mass of the nucleon and 
by S its polarization. The final (unpolarized) hadrons, with mass Mi, M2 and momenta Pi, 
P2, have invariant mass squared = (which we consider much smaller than the hard 
scale = —q^ > of the SIDIS process). The P^ = Pi+ P2 is the total momentum of the 
pair; we define also its relative momentum R = {Pi — P2)/2. The momentum transferred 
to the nucleon target is q = k — k'. The kinematics of the process is depicted in Fig. 1 (see 
also Ref. [26]). 

The component St of the target polarization is transverse to both the virtual-photon 
and target momenta q and P, respectively. Instead, Rt = 11— {R ■ Ph)Ph is orthogonal 
to Pfi but, up to subleading- twist corrections, it can be identified with its projection on the 
plane perpendicular to q and containing also St- The azimuthal angle (pn is the angle of 
Rt about the virtual-photon direction; similarly for the azimuthal angle (j)s of St- The 
explicit expressions are 

{qxk)- Rt {q x k) ■ {q x Rt) 

=77 1\ PI arccos — — — , 

\{q X k) ■ Rt\ \q x k\\q x Rt\ 

_ {qxk)-ST (g X fc) • (g X St) , . 

VS =y, , X I arccos — — — . (2.2) 

|(gxfc)-5T| |gxfc||gxST| 

We also define the polar angle d which is the angle between the direction of the back-to-back 
emission in the center-of-mass (cm) frame of the two hadrons, and the direction of P/^ in 
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Figure 1. Kinematics of the two-hadron semi-inclusive production. The azimuthal angles (j)^ 
of the dihadron, and (ps of the component St of the target-polarization, transverse to both the 
virtual-photon and target-nucleon momenta q and P, respectively, are evaluated in the virtual- 
photon-nucleon center-of-momentum frame. 



the photon-proton cm frame (see Fig. 1). We have 

1 



R\ = - ^Ml - 2{Ml + Ml) + {Ml - MIY/Ml , 
RT = Rsme. (2.3) 

As usual in SIDIS, we define also the following kinematic invariants 

Q2 p.q 2Mx 

P-^h _ ^ ^ ZX-Z2 , . 

z = s = zi + Z2 , C = , (2-5) 

P ■ q z 

where zi, Z2^ are the fractional energies carried by the two final hadrons. The invariant ^ 
can be shown to be a linear polynomial in cos0 [27]. 

To leading order in the couplings and leading twist, the differential cross section for 
the two-hadron SIDIS of an unpolarized lepton off a transversely polarized nucleon target 
contains only two nonvanishing structure functions: 

da 



dx dy dip dz d(pR dM'^ d cos 9 

^ Fuu + \ST\esin{c^R + <Ps)FS^^'^^+^A , (2.6) 



a2 ..2 




xyQ2 2(1 -e) 



where a is the fine structure constant and the structure functions F depend on x, Q^, z, 
cos9, and Mh- The first and second subscripts of F indicate the polarization of beam and 
target, respectively. Here, the target polarization refers to the virtual-photon direction; the 
conversion to the experimental polarization with respect to the lepton beam is straightfor- 
ward and given in Ref. [28]. The angle ip is the azimuthal angle of i' around the lepton beam 
axis with respect to an arbitrary fixed direction, which in case of a transversely polarized 
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target we choose to be the direction of S. The corresponding relation between tp and (ps is 
given in Ref. [28]; in DIS kinematics it turns out ~ d(ps- 

The ratio e of longitudinal and transverse photon flux in Eq. (2.6) is given by [29] 



1 - y - h^y'^ 



(2.7) 



so that the depolarization factors can be written as 

2(r^ = 1^ (1 - 2/ + + li'y") - (1 - 2/ + \y') - m , 

^(rr^^ = Y^(l-2/-hV) «(l-y) ^B{y). (2.8) 

The £ turns out to be also the ratio between the two depolarization factors: e = B{y) /A{y). 
Neglecting target-mass corrections, we will assume that in each experimental bin 

{A{y)) ^ A{{y)) , {B{y)) ^ B{{y)) , 
.^B{{y)) 



^ Cy . (2.9) 

^ <^ Q^, the structure functions c 

DiFFs [24, 27, 30]^ 



In the limit M? <^ Q^, the structure functions can be written as products of PDFs and 



Fuu = xJ2 el fUx; Q^) Dl{z, cose, Mh;Q^) , (2.10) 
g 

^jn(<AH+0s) =\^^^^elhl{x;Q')H<'^{z,cose,M,;Q') , (2.11) 

where Cq is the fractional charge of a parton with flavor q. The Df is the unpolarized DiFF 
describing the hadronization of a parton with flavor q into an unpolarized hadron pair 
plus anything else, averaging over the parton polarization. The H^'^ is a chiral-odd DiFF 
describing the correlation between the transverse polarization of the fragmenting parton 
with flavor q and the azimuthal orientation of the plane containing the momenta of the 
detected hadron pair. 

For Mfi ^ Q, the hadron pair can be assumed to be produced mainly in relative s or p 
waves, suggesting that the DiFFs can be conveniently expanded in partial waves. In the two- 
hadron cm frame, the relevant changes in the kinematics are summarized in Eq. (2.3). From 
the simple relation between and cos^, DiFFs can be expanded in Legendre polynomials 
in cos 9 as [27] 

Di ^i,ss+pp + -Di,sp cos e + D^ ppj (3 cos^ 0-1) , 



'4 

H<^H<,p + H<ppCos9. (2.12) 



^For some discussion of the case ~ , see Ref. [31] 
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All the cos ^—dependent terms disappear after integrating upon dcos9; they still vanish 
even if the 6 dependence of the acceptance is not complete but symmetric about 6 = Tr/2. 
Of the remaining terms, the subscript ss + pp refers to the unpolarized pair being created 
in a relative AL = state, while sp indicates the interference in \AL\ = 1. For simplicity, 
we will use the notation ss+pp = -^i since no ambiguity arises in the following. 

From the cross section of Eq. (2.6), and by inserting the structure functions of Eqs. (2.10), 
(2.11) with the approximation of Eq. (2.12), we obtain the following single-spin asymmetry 
(SSA) [24, 27, 32] 



n I ^^R ^'^^ ^ sin(0^ + (pg) {da^ — da^) ^ ^ J d cos 



,sin{</>^+0^) 



\St\ f dcj) d cos 9 {da^ + da^) 

^ \R\ elh'i{x;Q')Hl^^{z,Mj,;Q^) 

'Mj, Z,e^j1ix;Q^)Dliz,Mh;Q^) 

= ^SIDIS • 



/ dcosOFi 



uu 



(2.13) 



We are interested in the specific case of semi-inclusive production of 7r"'"7r pairs. Then, 
isospin symmetry and charge conjugation suggest [32, 33] 



Dl = Dl , Df = Df, Dl = Df 



TT< d T_r< U TT< d 

-Hi — -Ml — Hi , 



(2.14) 
(2.15) 



For a proton target, the SSA (2.13) simplifies to [33] 



\R 



Cy^H<^(z,Mh;Q^) K^ix;Q^) - \hf{x;Q^) 



X /r"(x; Q') D^z, Mh-, Q') + y^^ix; Q') Dl{z, M^; Q') 



+ \fi+\x)D{+\z,Mh;Q^ 



-1 



(2.16) 



and for a deuteron target to 



7 'irKspi^^ ^h-. g') k" {x- q2) + hf' [x- q2) 



Ag^iQ{x,y,z,Mh;Q'^) 
3 \R\ 

+ \f+\x-Q^)D{+\z,Mh;Q^)^ ' 



(2.17) 



where hf =h\-h\ and = + . 
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Eqs. (2.16) and (2.17) contains two sets of unknowns: the transversity hi (in various 
flavor combinations) and the DiFFs. Before the measurement by the Belle collaboration 
of the angular distribution of two pion pairs produced in e+e~ annihilations [34], the only 
information available on DiFFs were coming from model calculations in the context of the 
spectator approximation [24, 32, 35] (and, recently, also using the NJL-jet model [36]). The 
unpolarized Di was tuned to the Monte Carlo event generator [32] and the polarized H^^^ 
compared to the asymmetry measured by the HERMES collaboration in SIDIS on proton 
targets [37], the only available set of experimental data at that time [26]. 

The first analysis of the so-called Artru-Collins asymmetry [38] in e~^e~ annihilations 
by the Belle collaboration made possible a direct extraction of H^^^ for the production 
of tt^tt" pairs. In the absence of a measurement of the unpolarized e~^e~ cross section 
(planned at Belle in the near future), Di was parametrized to reproduce the two-hadron 
yield of the PYTHIA event generator, which is known to give a good description of data. 
Combining such a parametrization with the fit of the azimuthal asymmetry presented in 
Ref. [34], it was possible to extract for the first time the H^^^ [39]. 

The knowledge of DiFFs in Eq. (2.16) allowed us to get a glimpse of the combination 
_ directly from the HERMES data for Agims [33]. The effects produced by 

evolving DiFFs between the HERMES and Belle very different scales were properly included 
by using standard evolution equations in a collinear framework [25] and by implementing 
leading-order (LO) chiral-odd splitting functions in the HOPPET code [40]. Recently, the 
COMPASS collaboration has released new data for ^gimg on a proton target and for ^gj^jg 
on a deuteron target, with higher statistics and wider kinematic coverage [41]. Thus, the 
combination of SSA of Eqs. (2.16) and (2.17) makes it possible to separately parametrize 
each valence fiavor of the transversity distribution, which we present here for the first time 
(see also Ref. [42] for a first attempt to obtain a fiavor separation point by point). 

In the collinear framework, the dependence of the SSA on the momentum fraction x 
gets factorized from the dependence on {z,Mfi), as it is evident in Eqs. (2.16) and (2.17). 
This suggests that the dependence of the SSA on x comes only from the involved PDFs. 
Therefore, it is more convenient to study it by integrating the z- and M/^-dependence of 
DiFFs. Then, the actual combinations of transversity used in the analysis are, for the 
proton. 



X hl{x- Q2) = X h^^ (x; Q2) _ 1 xhf^ (x; Q^) 

_ ^SIDIs(^' 



n4Q')x/r"(^;Q') + inrf(Q^)x/f+'^(x;Q^) + i n,(g^) x/r'(x; Q^) 



es+s 



(2.18) 
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and, for the deuteron, 



x/if (x;Q2) = xh'^^{x;Q^) + xhf-{x;Q^) 

_ ^SIDIs(^' Q^) ^ 

CynliQ^) 3 

where 



2, , nd(Q2)\ ns{Q^) 



+ 



xf^+%x;Q-') 













■^min min 




r^max r-^^h. max 


n\{Q') = 






-^min -^'-^h min 



dzdMhDl{z,Mh;Q^) , 
dzdMj-^^H^^l{z,Mh;Q') 



(2.19) 

(2.20) 
(2.21) 



Using the unpolarized PDFs from the MSTW08LO set [43], early explorations were pre- 
sented in Ref. [44]. 



3 Fitting procedure 

Here, we describe our fitting procedure to obtain the valence transversity distribution func- 
tions for up and down quarks. We discuss first the choice of the functional form. 
The main theoretical constraint we have is Soffer's inequality [45] (see also [46]) 

2\hlix;Q^)\ < \fl{x;Q^) + glix;Q^)\^2SB''ix;Q'') . (3.1) 

We impose this condition by multiplying the functional form by the corresponding Soffer 
bound at the starting scale of the parameterization. If the Soffer bound is fulfilled at some 
initial Qq, it will hold also at higher > Qq [47, 48]. The implementation of the Soffer 
bound depends on the choice of the unpolarized and helicity PDFs. We use the MSTW08 
set [43] for the unpolarized PDF, combined to the DSSV parameterization [7] for the helicity 
distribution, at the scale of Qg = 1 GeV^. Our analysis was carried out at LO in as- To be 
as consistent as possible, we decided to use the MSTW08LO set for /i and the DSSV set 
for gi, even if the DSSV fit provides only a NLO parametrization of gi. For convenience, 
in App. A we list the explicit form of SB'^. The result for the Soffer bound is affected by 
an error coming mainly from the uncertainty in the knowledge of the helicity PDF gi. We 
checked that at the explored hard scales this error is much smaller than the experimental 
errors on ^siDls data and the statistical error on the DiFF parametrization; hence, we will 
neglect it. 

The Soffer bound is valid for each quark and antiquark. Since we need to parametrize 
the transversity valence combinations for up and down quarks, we have necessarily to con- 
strain it by taking the sum of Soffer bounds for both quarks and antiquarks. This likely 
leads to a loose bound, especially at low x. In particular, due to the divergent behavior of 
PDFs the "valence" Soffer bound is not even integrable in the range x G [0, 1], which would 
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result in a divergent tensor charge. Given the chosen analytical form of the PDFs at low x, 
the SB^ has to be multiphed by at least x°-^^^'^^. 

Based on the above considerations, we adopted the following functional form for the 
valence transversity distributions at = 1 GeV^: 

xhf'{x;Ql) =tanh x^^"^ (Ag+Bgx + Cqx'^+Dgx^) x SB«(x; Qg) +3; SB«(x; Qg) . (3.2) 

The hyperbolic tangent is such that the Soffer bound is always fulfilled. The functional 
form is very flexible and can contain up to three nodes. The low- a; behavior is however 
determined by the x^/^ term, which is imposed by hand. Present fixed-target data do not 
allow to constrain it. 

The fit, and in particular the error analysis, was carried out in two different ways: using 
the standard Hessian method and using a Monte Carlo approach. As usual, we remind the 
reader that both methods are suitable to estimate the errors of statistical nature only, 
assuming a specific choice of the theoretical function. Care must be taken especially when 
using error bands outside the region where data exist. 

The standard fitting procedure consists in minimizing the usual function, defined as 

(xi hl^ZSxi; Q2) _ X, h^/Zoixi, Ql, {p})) ^ 
X\{P}) = E ^ 7 JJ, ^2 ' (3-3) 



where the sum runs over the experimental points, the expressions for x hF^^ are listed in 
Eqs. (2.18) and (2.19), and {p} denotes the vector of parameters. The evolution of the 
functional form (3.2) to the values Qf for each data bin has been implemented using the 
HOPPET code [40], set to the MS renormalization scheme and modified to include also 
the LO chiral-odd splitting functions needed for transversity evolution. The input value 
for the running coupling constant at Qq = 1 GeV^ is chosen to be the best-fit value of 
the MSTW08LO set, i.e. df{Ql) = 0.13939. It is true for all the evolved quantities of 
our analysis, including the Soffer bound. The minimization has been carried out using the 
MINUIT code and led to a vector of best-fit parameters, {po} (and a covariance matrix). 

The standard method allows to compute the errors on any theoretical quantity under 
the assumption that the parameter dependence of can be approximated by a quadratic 
expansion around the minimum, and the parameter dependence of the theoretical quantity 
can be approximated by a linear expansion around the minimum. 

For the standard method, the error on the extracted transversity was estimated using 
the formula 

A/ii(x,(3 ■,{p])) = > ^ Covij , (3.4) 

/ ^ dpi {po} dpj {po) 

where Np is the number of parameters. The covariance matrix Covjj has been obtained 
using the condition Ax^ = 1. Therefore, within the limits of applicability of the standard 
approach, the obtained error band corresponds to the lo" or 68% confidence level. In typical 
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PDF global fits, often the value of is increased of one or even two orders of magnitude, 
with a corresponding increase in the error estimate (see, e.g., [7]). In the present analysis, 
we find no need of such an increase, as demonstrated by the agreement with the Monte 
Carlo approach. 

The Monte Carlo approach does not rely on the assumptions of a quadratic dependence 
of and a linear expansion of the theoretical quantity around {po}, respectively. In 
our case, the need of such an approach is essential whenever the minimization pushes the 
theoretical functions towards their upper or lower bounds, where it is not possible to assume 
a simple linear expansion in the parameters. 

For the implementation of this approach, we took inspiration from the work of the 
NNPDF collaboration (see, e.g., [49-51]), although our results are not based on a neural- 
network fit. The approach consists in creating N replicas of the data points. In each replica 
(denoted by the index r), each data point i is shifted by a Gaussian noise with the same 
variance as the measurement. Each replica, therefore, represents a possible outcome of an 
independent experimental measurement, which we denote by h^^^ [xi] Q?). The number of 
replicas is chosen so that the mean and standard deviation of the set of replicas accurately 
reproduces the original data points. In our case, we have found that 100 replicas are 
sufficient. 

The standard minimization procedure is applied to each replica separately, by mini- 
mizing the following error function^ 



resulting in N different vectors of best-fit parameter values, {por}; f = 1, • . . iV. These 
parameter vectors can be used to produce values for any theoretical quantity. The 
theoretical outcomes can have any distribution, not necessarily Gaussian. For non-Gaussian 
distributions, the la confidence interval is in general different from the 68% interval. Both 
of them can be easily computed from the N theoretical outcomes. For instance, for the 
68% interval we simply take for each experimental point i the A'^ values and we reject the 
largest and the lowest 16% of them. 

Although the minimization is performed on the function defined in Eq. (3.5), the agree- 
ment of the A^ theoretical outcomes with the original data is better expressed in terms of 
the original function defined in Eq. (3.3), i.e. with respect to the original data set with- 
out the Gaussian noise. If the model is able to give a good description of the data, the 
distribution of the A'^ values of x^/d.o.f. should be peaked at around one. In real situations, 
the rigidity of the model shifts the position of the peak to higher values of x^/d.o.f.. 

In our case, we determined A^ = 100 best-fit parameter vectors and we used them to 
produce 100 curves for the up and down valence transversity. Each one of the resulting 

^Note that the error for each rephca is taken to be equal to the error on the original data points. This 
is consistent with the fact that the variance of the replicas should reproduce the variance of the original 
data points. 
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HERMES 


data 






X 


V 


g2[GeV2] 


^SIDIS 


/ii"" - /if" /4 


0.033 


0.734 


1.232 


0.015 ±0.010 


0.086 ±0.061 


0.047 


0.659 


1.604 


0.002 ±0.011 


0.010 ±0.054 


0.068 


0.630 


2.214 


0.035 ±0.011 


0.167 ±0.069 


0.133 


0.592 


4.031 


0.020 ±0.010 


0.092 ±0.054 




COMPASS 


proton 


data 




X 




O^fCeV^l 


r\ O TT^ TO 


/i"" - /if" /4 


\J .\j\J\jO 




1.232 


n nO(^ _L n non 




01 0^ 

\J.\J ±\jO 






u.u±u zc u.u±u 


U.UOO ZC \J.\JOiJ 


0.0164 




1.744 


0.015 ±0.013 


0.057 ± 0.049 


0.1330 




2.094 


0.008 ±0.010 


0.031 ±0.039 


0.0398 




2.802 


0.027 ±0.011 


0.107 ±0.049 


0.0626 




4.342 


0.029 ±0.014 


0.118 ±0.060 


0.1006 




6.854 


0.051 ±0.016 


0.208 ± 0.079 


0.1613 




10.72 


0.108 ±0.023 


0.42 ±0.12 


0.2801 




21.98 


0.080 ± 0.033 


0.24 ±0.11 




COMPASS 


deuteroii 


data 




X 




O^fCeV^l 

[VJ V J 


A\ O TT^ TO 


, , , 


0.0064 




1.253 


\j _i_ 








-L . (J WO 




U.Wt: —I— U.J-Zj 


0.0163 




1.792 


0.028 ±0.010 


0.28 ±0.11 


0.0253 




2.266 


-0.005 ±0.009 


-0.051 ±0.094 


0.0396 




3.350 


0.006 ±0.011 


0.06 ±0.12 


0.0623 




5.406 


-0.006 ±0.014 


-0.06 ±0.14 


0.0996 




8.890 


-0.029 ±0.019 


-0.30 ±0.20 


0.1597 




15.65 


-0.017 ±0.030 


-0.16 ±0.28 


0.2801 




33.22 


0.078 ±0.054 


0.50 ±0.36 



Table 1. HERMES data for tt+tt^ production in SIDIS off a transversly polarized proton [26] 
and COMPASS data for the same process off a transversly polarized proton and deuteron [41]. The 
last column shows the combinations of valence transversities obtained using Eqs. (2.18) and (2.19). 

curves respects the Soffer bound by construction. The results are discussed in the next 
section. 

4 Results and discussion 

In the following, we discuss the results obtained by fitting the expressions of Eqs. (2.18) 
and (2.19) when inserting the HERMES and COMPASS measurements for the single-spin 
asymmetries ^gi^jg and ^gjoig on vr+vr" SIDIS production off transversely polarized pro- 
ton and deuteron targets, respectively. By combining the two fits, we can determine for 
each valence flavor and the vector of fitting parameters that gives the corresponding 
transversity distribution, according to Eq. (3.2). 

As discussed in the previous section, the error analysis has been performed in two ways: 
using the standard Hessian method summarized in Eqs. (3.3) and (3.4), or the Monte Carlo 
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approach by fitting = 100 replicas of the experimental points according to Eq. (3.5). 
For each strategy, we explored three different scenarios in the parametrization (3.2) of the 
transversity distribution: 

• the rigid scenario, described by the choice Cu = = Du = = 0, i.e. with only 4 
free parameters; 

• the flexible scenario, with = = (6 free parameters); 

• the extra-flexible scenario, with all 8 free parameters. 

4.1 Experimental data 

In Tab. 1, we list the data for the asymmetries ^gj£)jg and ^gj^jg of Eqs. (2.16) and (2.17), 
as they were measured by the HERMES [26] and COMPASS [41] collaborations in the 'k~^tt~ 
SIDIS production off transversely polarized proton and deuteron targets, respectively. The 
first three columns indicate the average values of the corresponding kinematic variables in 
each experimental bin. The indicated errors include statistical and systematic contributions 
added in quadrature. The depolarization factor Cy in the expression of the asymmetries 
depends on the average y according to Eq. (2.9). In its analysis, the COMPASS collabora- 
tion already divided the Cy factor out of the measured cross section; as such, the SSA does 
no longer depend on y and, correspondingly, there are no y values in the second column 
of Tab. 1 for COMPASS. Consistently, we have used Cy = 1 when fitting the COMPASS 
data. 

The last column in Tab. 1 contains the values of the combination in Eq. (2.18) for 
the proton target, and of Eq. (2.19) for the deuteron target, when the corresponding ex- 
perimental values for the SSA are inserted in ^gj^jg and ^gj^jg, respectively. As already 
anticipated in the previous sections, for the unpolarized PDFs we adopted the MSTW08LO 
set [43]. The remaining ingredients in Eqs. (2.18) and (2.19) are the Uq and defined in 
Eqs. (2.20) and (2.21), where the DiFFs D\ and -fffj^ are parametrized as in Ref. [39]. The 
integrals are evaluated according to the appropriate experimental cuts: 0.2 < z < 1 and 
0.5 GeV < Mh<l GeV for HERMES, 0.2 < z < 1 and 0.29 GeV < Mh < 1.29 GeV for 
COMPASS. In Tab. 2, the results are given for the relevant flavors at the average scales 
for each experimental bin i. The statistical error is indicated only for n^, since the large 
statistics achievable in the Monte Carlo simulation of the unpolarized e~^e~ — )• (7r''"7r^)X 
cross section makes the error of Uq negligible. 

In the standard Hessian method, the best fit parameters and their \a error (corre- 
sponding to = 1) at the initial scale Qq = 1 GeV^ are given in Tab. 3. The x^/d.o.f. 
is 1.23 for the rigid scenario, 1.12 for the flexible scenario, and 1.26 for the extra-flexible 
scenario. For the 100 replicas, the average x^/d.o.f. are 1.35, 1.56 and 1.86, respectively. 

4.2 Fitting results 

In Fig. 2, the points with error bars represent the transversity combinations for proton (left 
panel) and deuteron target (right panel) quoted in the last column of Tab. 1. The central 
value of our best-fit result in the standard approach with the flexible scenario is given by the 
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Table 2. The integrated DiFFs according to Eqs. (2.20) and (2.21). The error has been computed 
at the average for each indicated experimental bin. 

thick solid line, and it is in good agreement with the data. The error band is the outcome 
of the merging of all the straight lines connecting the statistical error bars of the fit for each 
experimental point. The other scenarios do not show significant qualitative differences in 
the range where data exist. 

In Fig. 3, we show the same comparison in the same conditions as in the previous figure, 
but for the Monte Carlo approach. The band now represents the result of the 68% of all 
replicas, obtained by rejecting the largest 16% and the lowest 16% of the replicas' values in 
each X point. We observe no substantial difference between the standard and Monte Carlo 
approaches. 

The resulting transversity distribution is plotted in Fig. 4. The left panel displays the 
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Table 3. Best-fit parameters and values obtained in the standard approach for the three 
scenarios described in the text and based on Eq. (3.2). 



X h,"v(x)-x hi^vCxjM X h,"v(x)+x hi'^v(x) 




Figure 2. Tire combinations of Eq. (2.18), left panel, and Eq. (2.19), right panel. The squares 
and triangles are obtained from the COMPASS and HERMES data, respectively (the values are 
indicated in the last column of Tab. 1). The thick solid line indicates the central value of the 
best-fit result in the standard approach with the fl.exible scenario (see text). The error band is the 
outcome of the merging of all the straight lines connecting the statistical error bars of the fit for 
each experimental point. 

q = Uv contribution in Eq. (3.2), while q = is in the right one. From top to bottom row, 
the results for the rigid, flexible, and extra-flexible, scenarios are shown, respectively. For 
each panel, the outcome in the standard approach with the Hessian method is represented 
by the uncertainty band with solid boundaries, the central thick solid line visualizing the 
central value. The partially overlapping band with dashed boundaries is the outcome when 
adopting the Monte Carlo approach, where the band width corresponds to the 68% of all 
the 100 replicas, again produced as before by rejecting the largest 16% and the lowest 
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Figure 3. Same observables and data symbols as in the previous figure. The uncertainty band 
represents in the Monte Carlo approach the selected 68% of all fitting replicas (see text). 



16% among the replicas' values in each x point. As such, the set of selected replicas in 
the 68% band can change in each different x point; consequently, the band itself can show 
some irregular wiggles. For sake of comparison, each panel displays also the corresponding 
results for the only other existing parametrization available [15], depicted as a band with 
short-dashed boundaries. Since the latter was extracted at the scale = 2.4 GeV^, our 
results are properly evolved at the same scale. Finally, the dark thick solid lines indicate 
the Soffer bound, also evolved at the same scale = 2.4 GeV^ (using LO evolution as in 
the rest of the analysis). 

For the flexible scenario (middle row of Fig. 4), the uncertainty bands in the standard 
and Monte Carlo approaches are quite similar. The main difference is that in the former 
case the boundaries of the band can occasionally cross the Soffer bound. This is due to the 
fact that the assumed quadratic dependence of the parameters around its minimum 

is not a reliable one, when getting close to the bounds. On the contrary, in the Monte Carlo 
approach each replica is built such that it never violates the Soffer bound; the resulting 
68% band is always within those limits. 

For the valence up contribution (left panel), the standard approach tends to saturate 
the Soffer bound at x ~ 0.4 (outside the range where data exist). In the Monte Carlo 
approach, some of the replicas saturate the bound already at lower values. However, there 
are also a few replicas that do not saturate the bound at all, or even saturate the lower 
Soffer bound. These replicas typically fall outside the 68% band drawn in the figures. 
Nevertheless, they can still have a good when compared to the data. 

For the down valence contribution (right panel), both approaches saturate the lower 
limit of the Soffer bound already at x ~ 0.1, i.e. in a region where data exist. This 
behavior is driven by the data, in particular by the bins number 7 and 8 in the deuteron 
measurement. No such trend is evident in the corresponding single-hadron measurement 
of the Collins effect, from which the other parametrization of Ref. [15] is extracted. As a 
matter of fact, this is the only source of significant discrepancy between the two extractions, 
which otherwise show a high level of compatibility despite the fact that they are obtained 
from very different procedures. Note that if the Soffer bound is saturated at some scale, it 
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Figure 4. The up (left) and down (right) valence transversities coming from the present analysis 
evolved to = 2.4 GeV^. From top row to bottom, results with the rigid, flexible, and extra-flexible 
scenarios are shown, respectively. The dark thick solid lines are the SofFer bound. The uncertainty 
band with solid boundaries is the best fit in the standard approach at la, whose central value is 
given by the central thick solid line. The uncertainty band with dashed boundaries is the 68% of 
all fitting replicas obtained in the Monte Carlo approach. As a comparison, the uncertainty band 
with short-dashed boundaries is the transversity extraction from the Collins effect [15]. 



is likely to be significantly violated at a lower scale [46]. Therefore, if we want to maintain 
the validity of the Soffer bound at < 1 GeV^, we would expect transversity to be clearly 
below the Soffer bound at > 1 GeV^. In fact, in our analysis with the Monte Carlo 
approach there are a few replicas that do not saturate the Soffer bound. They fall outside 
the 68% band drawn in the figure, but they are still compatible with the data due to the 
large experimental error bars (this is true in particular for the deuteron bins number 7 and 
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8). Therefore, at present we cannot conclude that the the Soffer bound is saturated or 
violated, even though the fit seems to point in that direction. We mention that interesting 
speculations concerning violations of the Soffer bound were presented in Ref. [52]. 

At low X, the functional form in the flexible scenario tends to zero by construction, and 
similarly in all other scenarios. However, the behavior down to x ~ 0.005 is driven by data. 
In fact, the functional form can have up to two nodes in x G [0,1]. In the Monte Carlo 
approach, most of the replicas for x/i^" have no node, while for x/i^" have one or even two 
nodes. 

In the rigid scenario (upper row in Fig. 4), most of the features are similar to the flexible 
scenario in the region where data exist, but there are some differences outside that range. 
For the valence up quark (left panel), both standard and Monte Carlo approaches give 
uncertainty bands that saturate the Soffer bound at higher x, almost completely overlapping 
with the result obtained from the Collins effect. For the valence down quark (right panel), 
the trend is very similar to the flexible scenario: the parametrization saturates the lower 
Soffer bound at x > 0.1. This demonstrates that this unexpected behavior is not an artefact 
of the functional form, but is due to the experimental data. In the Monte Carlo approach, 
the majority of replicas show the same behavior, but a few ones (falling outside the 68% 
band) do not, again as in the flexible scenario. 

Finally, for the extra-flexible scenario (bottom row in Fig. 4) the distinction between 
regions with and without experimental data is even more clear. Where there are data, 
the results are highly compatible with the other scenarios. But at a; > 0.4 for the 
case (left panel), the uncertainty band in the standard approach substantially violates the 
Soffer bound. This is an artefact of the assumptions used in error propagation, together 
with the lack of data at high x. In the Monte Carlo approach, the replicas entirely fill 
the area between the upper and lower Soffer bound, both for the up quark at x > 0.4 and 
for the down quark at x > 0.25. This is an explicit visualization of the realistic degree of 
uncertainty about transversity in the x range where there are no experimental data points. 

We have explored other scenarios for the transversity functional form. We tried different 
arguments in the hyperbolic tangent, specifically different powers in the first factor, but with 
no significant change for x > 0.01. Using the x^/^ factor, the error band at low x becomes 
considerably wider. The valence transversities remain integrable, but stable values of the 
tensor charge can be reached only by pushing the lower limit in the integral to extremely 
small values x < 10^^*^ (see below). 

A qualitative comparison between the results of the present work and the available 
model predictions can be done using the results collected in Ref. [53]. In particular, we 
note that the extracted transversity for the up quark is smaller than most of the model 
calculations at intermediate x S [0.1, 0.2], while it is larger at lower x (x ~ 0.01). The down 
transversity is much larger in absolute value than all model calculations at intermediate x 
(as observed before, this is due to the deuterium data points), while the error band is too 
large to draw any conclusion at lower x. 

Transversity is directly related to the tensor charge, a fundamental quantity of hadrons 
at the same level as the vector, axial, and scalar charges. The tensor charge remains at the 
moment largely unconstrained. It can be directly compared with lattice QCD predictions 
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-0 26 ± 17 
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MC rigid 


0.30 ±0.07 


-0.22 ±0.11 


0.56 ±0.12 


-0.08 ±0.27 
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0.29 ±0.13 


-0.26 ±0.22 


0.72 ±0.24 


-0.33 ±0.61 


MC flex. 


0.32 ±0.09 


-0.24 ±0.11 


0.77 ±0.22 


-0.45 ±0.48 
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0.32 ±0.12 


-0.25 ±0.15 


0.61 ±0.40 


-0.16 ±0.44 
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0.34 ±0.10 


-0.20 ±0.14 


0.68 ±0.22 


-0.12 ±0.69 



Table 4. Table of the results for the tensor charge at Qq = 1 GeV, truncated in the range where data 
exist (second and third column) and extended to the whole x range (third and fourth column) . The 
results are given for the standard and Monte Carlo approach and for the three scenarios considered 
in the fit. 



(see, e.g., Refs. [54, 55]) or models (see, e.g., Refs. [56-60]). There is no sum rule related 
to the tensor current, due to the property of the anomalous dimensions governing the QCD 
evolution of transversity. The contribution of a flavor q to the tensor charge is defined as 

HQ^) = J dxh\^{x;Q'^) . (4.1) 

The region of validity of our fit is restricted to the experimental data range. We 
can therefore give a reliable estimate for the tensor charge truncated to the interval x G 
[0.0064,0.28]. In the first two columns from left of Tab. 4 we list the results obtained in 
the different approaches and scenarios. We tried also to extend the range of integration 
outside the experimental data to x G [0, 1]. The result is heavily influenced by the adopted 
functional form, in particular by the low-x exponent. Nevertheless, we quote our result in 
the last two columns of Tab. 4. 

Our results for the tensor charges in the flexible scenario are slightly larger in absolute 
value compared to the ones in Ref. [15]. They are compatible within errors. Results 
obtained in several models can be found in, e.g., Refs. [15, 61, 62]. As emphasized in 
Ref. [62], care must be taken when comparing results at different scales. The ratio between 
the up and down tensor charges is scale invariant. The values we obtain in all our scenarios 
are compatible with all models within the large errors. In order to better determine the 
tensor charge, more data at high and low x are needed. 

5 Conclusions and Outlook 

The transversity parton distribution function (PDF) is an essential missing piece of our 
knowledge on the proton at leading- twist. It is a chiral-odd object, whose dependence 
obeys the non-singlet QCD evolution. Its integral over x is related to the nucleon tensor 
charge. Positivity bounds constrain its absolute value to be smaller than the absolute value 
of the number density and helicity, the so-called Soffer bound. Due to its chiral-odd nature, 
transversity cannot be accessed in fully inclusive deep-inelastic scattering. It is however 
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possible to access it in two-particle-inclusive DIS [22-24] in combination with Diliadron 
Fragmentation Functions (DiFFs). 

In this paper, we have obtained for the first time the parameterization of the up and 
down valence transversities based on a collinear framework, using data for vt+tt" semi- 
inclusive DIS off transversely polarized targets from the HERMES and COMPASS collab- 
orations [26, 41], combined with the Belle data on almost back-to-back emission of two 
7r+7r~ pairs in e~^e~ annihilations [34]. We have explored different scenarios for the func- 
tional form, all subject to the theoretical constraint of the Soffer bound [45, 46]. We have 
also performed the error analysis in two independent ways. The first one is a standard 
one based on the Hessian method. The second one is based on the random generation 
of a large number of replicas of the experimental points, and on the fit of each of these 
replicas, producing an envelope of trajectories whose spread is the generalization of the la 
uncertainty band when the distribution is not necessarily a Gaussian. As such, the second 
method is more reliable particularly when the fitting curves hit the Soffer bound, and the 
function cannot be expected to have the quadratic dependence on the fit parameters as 
required by the Hessian method. Nevertheless, in the kinematical range of the experimental 
measurements the two methods give almost overlapping results in all explored scenarios. 

In the range where data exist, our results are compatible with the only other existing 
parametrization of transversity, which is determined from the Collins effect in single-hadron 
SIDIS off transversely polarized targets [15]. The only source of discrepancy lies in the range 
0-1 i$ ^ i$ 0-16 for the valence down quark, where two experimental data for the deuteron 
target drive our fitting curves to saturate the lower Soffer bound. However, the large error 
bars of these two points prevent us from drawing any conclusion about a possible violation 
of the Soffer inequality. 

Outside the kinematical range of experiments, the lack of data reflects itself in a large 
uncertainty in the parametrization. For the extra-flexible scenario and the error analysis 
based on the random approach, the replicas take all the available values between the upper 
and lower Soffer bound at large x. This illustrates in a very effective way the need for new 
large— rr data in order to reduce the degree of uncertainty in the knowledge of transversity. 

In the near future, more data are expected from the HERMES and COMPASS collab- 
orations. They will include also different types of hadron pairs (e.g., Ktt), which should 
allow us to improve the flavour separation of transversity. Two-particle inclusive DIS will be 
measured also at JLab in the future, which should considerably increase our knowledge of 
transversity at high x. Finally, invaluable information will come also from polarized proton- 
proton collisions [63]: data are expected from the PHENIX and STAR collaborations (see, 
e.g., [64]). 
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A Appendix 



For convenience, we reproduce here the exphcit analytic forms of the Soffer bound used in 
our analysis, which is based on the MSTW08LO set [43] for the unpolarized PDF, combined 
to the DSSV parameterization [7] for the helicity distribution. The equations hold at Qq = 1 
GeV2. 
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